Quantum spin Hall effect and spin-charge separation in kagome lattice 
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A two-dimensional kagome lattice has been theoretically investigated within a simple tight-binding 
model, which includes the nearest neighbor hopping term and the intrinsic spin-orbit interaction 
between the next nearest neighbors. By using the topological winding properties of the spin-edge 
states on the complex-energy Riemann surface, the spin Hall conductance is obtained to be quantized 
as —e/2n {e/2n) in insulating phases. This result keeps consistent with the numerical linear-response 
calculation and the Z2 topological invariance analysis. When the sample boundaries are connected 
in twist, by which two defects with vr flux are introduced, we obtain the spin-charge separated 
solitons at 1/3 (or 2/3) filling. 
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Over the last two decades the topological band insu- 
lators (TBIs) have been a subject of great interest in 
condensed matter fieldii^. Different from the normal 
band insulators, the TBIs have a prominent characteris- 
tic, which is the necessary presence of gapless edge states 
on the sample boundarjs^ii. An early TBI model was 
proposed by Haldane^. Therein it was shown that the 
gapless edge states result in a remarkable feature of TBIs 
by showing a quantum Hall effect in absence of a macro- 
scopic magnetic field. Besides the Haldane model, several 
other lattice models have been proposed to be the quan- 
tum Hall TBIs, which include the two-dimensional (2D) 
spin-chiral kagome lattice^ii and the three-dimensional 
(3D) distorted fee lattice^i^. All the quantum Hall TBIs 
rely on the breaking of time-reversal symmetry (TRS) . 

Recently, Kane and Mele generalized the spinless Hal- 
dane model to a spin one by adding an intrinsic spin- 
orbit interaction (SOI)i£iii. The TRS is conserved in 
the Kane-Mele model and the gapless spin edge states in 
this model result in a quantum spin Hall effect. These 
TBIs like the Kane-Mele model keeping TRS are different 
from the quantum Hall TBIs, and we call them the quan- 
tum spin Hall TBIs. The quantum spin Hall TBIs have 
received considerable attention. One impressive example 
is that only one year later after its theoretical prediction 
to be a quantum spin Hall TBI in 2006^^, it was proved 
in experiment that HgTe is an actual onei^. 

One special character of the quantum spin Hall TBIs 
was recently attributed to their spin-charge separated 
excitationsi^ii^ in the presence of a tt flux. This attri- 
bution is motivated by the recent advance in studying 
2D fractionalized quasiparticle a^^'^^ , and is a straight- 
forward result when, like what Kane and Mele^^ have 
dealt with Haldane's TBI, considering spins of the 2D 
edge soliton. The separate spinon, holon, and chargeon 
obey Bose statistics and the experimental measurement 
of these soliton excitations would provide an undoubted 
verification of the Z2 topological properties of the quan- 
tum spin Hall TBIs. At present, besides the necessity 
for further studies to gain more insights into the nature 




FIG. 1: (Color online). Schematic picture of the 2D kagome 
lattice. The dashed lines represent the Wigner-Seitz unit cell, 
which contains three independent sites (A, B, C). 



of spin-charge separation and its connection to the other 
topological phenomena, obviously, it is also important 
to identify and study various model systems that exhibit 
the phenomenon of spin-charge separation. Motivated by 
this observation, as well as by the recent attention on the 
layered metal oxides as possible candidates for the quan- 
tum spin Hall TBI^'*, in this paper, we study the quantum 
spin Hall effect and construct spin-charge separated edge 
solitons in 2D kagome lattice. Different from the previ- 
ously studied kagome TBl'^''^'ia whearein the presence of 
ferromagnetic spin chirality breaks TRS, here TRS per- 
sists and the quantum spin Hall effect occurs due to the 
intrinsic SOI. By using the bulk linear-response theory, 
as well as the topological winding numbers of the spin- 
edge states on the complex-energy Riemann surface, we 
obtain the spin Hall conductance (SHC) a^y. The quan- 
tized value of a^y is — e/27r (e/27r) at 1/3 (2/3) fiUing. 
Then, we construct spin-charge separated edge solitons 
by introducing tt fluxes with a method similar to that in 
Ref.^^. The quantum statistics of these solitons is also 
discussed. 

Consider the tight-binding model for independent elec- 
trons on the 2D kagome lattice (Fig. [T]). The spin- 



independent part of the Hamiltonian is given by 
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where tij —t 
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the hopping ampUtude between the 
nearest- neighbor hnk (ij) and {cic) is the creation 
(annihilation) operator of an electron with spin a (up or 
down) on lattice site i. For simplicity, we choose t——l as 
the energy unit and the distance a between the nearest 
sites as the length unit throughout this paper. 

The Hamiltonian ^ can be diagonalized in the mo- 
mentum space as 



Wo = 5^V',^<,ffo(k)Vk., 
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where the electron field operator 

'0kCT=(cAk(T, c_Bk(T, cckcr)"^ includes the three lattice 
sites (A, B,C) in the Wigner-Seitz unit cell shown in 
Fig. 1. i?o(k) is a 3x3 spinless matrix given by 

(0 2cos(k-ai) 2cos(k-a3) 

2cos(k-ai) 2cos(k-a2) 

2cos(k-a3) 2cos(k-a2) 

■(3) 

where ai = (-l/2, - V3/2), a2 = (l,0), and 
a3=(— 1/2, -\/3/2) represent the displacements in a 
unit cell from A to B site, from B to C site, and from 
C to A site, respectively. In this notation, the first 
Brillouin zone (BZ) is a hexagon with the corners of 
K=± (2^/3) ai, ± (2^/3) as, ± (2^/3) aa. 

The energy spectrum for spinless Hamiltonian i?o(k) 

is characterized by one dispersionless flat band (e^^''=2), 
which reflects the fact that the 2D kagome lattice 
is a line graph of the honeycomb structure^^, and 



two dispersive bands, e 



(0) 
2(3)k 



= -l=F\/4&k - 3 with 



6k=X]i=i '^o^^ (k • a^). These two dispersive bands touch 
at the corners (K-points) of the BZ and exhibit Dirac- 

type energy spectra, e^^^^^^—{—l^\/3\k—K.\), which im- 
plies a "particle-hole" symmetry with respect to the 
Fermi energy eF= — 1- The corresponding eigenstates 
of -ffo(k) are given by 



G„k (9lk, '?2k, 93k) 



(4) 



where the expressions of the components qtk and the nor- 
malized factor G„(k) for each band are given in Table I. 



TABLE I: The expressions for the coefficients in Eq. Q with 
Xi=\^•^.i . 



qik 
q2k 

G 



1 r (0)2 



4 cos X2] 
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nk 



(0)2 



^nfc ''O^ X\ +2 COS X2 COS X3 
^nk COS X3 +2 COS X2 COS Xl 

[4bfe-36if ]cos^:c2+6(6fe-l)6i°,> 



FIG. 2: (Color online), (a) Energy spectrum (solid curves) of 
the 2D kagome lattice along the high-symmetry lines in the 
BZ with intrinsic spin-orbit couplings Aso^O.lt. There are 
two band gaps appearing with gap width A and S. For com- 
parison, we also draw the energy spectrum without intrinsic 
spin-orbit couplings (dotted lines), (b) The corresponding BZ 
of the 2D kagome lattice, (c) The SHC a^y as a function of 
the Fermi energy ep- The shaded areas correspond to the 
bulk gaps. 



Then, we introduce the intrinsic SOI term, which, ac- 
cording to the symmetry of the kagome lattice, takes the 
formia^ 



^so 
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so 
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Here Aso represents the SOI strength, s is the vector of 
Pauli spin matrices, i and j are next-nearest neighbors, 
and d}. and d?- are the vectors along the two bonds that 
connect i to j. Taking the Fourier transform, we have 



kcr 



V'kCT^so(k)V'kCT 



with 



Hso(k) -±2A, 



so 



( ° 

—i cos(k-bi) 
y icos(k-b3) 



icos(k-bi) 


—i cos(k-b2) 



=a2 



-i cos(k-b3) 
i cos(k-b2) 

(6) 

ai, and the -l-(— ) 



where bi—a3 —a2, b2=ai— 33, b3 
sign refers to spin up (down) electrons. 

Inclusion of the intrinsic SOI in the Hamiltonian makes 
the appearance of the eigenstates |M„k) and eigenener- 
gics e„k very tedious except for at some high-symmetry 
k points. Instead of writing their explicit forms, here 
we show in Fig. Ufa) (solid curves) the numerically 
calculated energy spectrum for the total Hamiltonian 
H—Ti-o+Hso along the high-symmetry lines (F^K, 
K^M, and M^F) in the BZ. The SOI coefRcient is cho- 
sen to be Aso=0.1i. For comparison we also plot in Fig. 



[SJa) (dashed curves) the energy spectrum in the absence 
of SOI (Aso=0). One can see that while the spin degen- 
eracies is not hfted by the presence of the intrinsic SOI, 
nevertheless, the Dirac contacts of the lower and mid- 
dle bands at the K point are removed and there evolves 
a direct gap between these two bands. The amplitude 
of this gap turns out to be A=4\/3|Aso|- Similarly, the 
original contact at the T point between the middle and 
upper (flat) bands is also lifted by the presence of SOI 
and a gap with amplitude 6 is opened. However, this gap 
is an indirect one, i.e., the middle-band maximum and 
upper-band minimum are not at the same fc-point. 

To see the behavior of the system in the insulating 
state, we have calculated the SHC using the following 
Kubo formula3^ 

= E [/(£nk)-/(e„'k)] (7) 

Im(u„k|^ {vx,Sz} |u«'k)(u n'kl'^y I'^nk) 
(enk-e„'k) +r] 

where v (k)=dH(k)/ndk and H(k)=Ho{'k) + ffso(k). 
The calculated result at zero temperature is shown in 
Fig. Etc) by varying the Fermi energy ep with the SOI 
coefficient Aso=0.1i. From Fig. 2(c), one can see that 
initially the SHC a^y decreases as the filling factors of 
the (spin-degenerate) lower band increases, arriving at 
the minimum value — (e/27r) at e^=— 1.35 (=— 1— A/2), 
a value corresponding to the top of the lower band. Then, 
as the Fermi energy ep continues to vary in the first 
gap region (shaded area wherein — 1.35^6;-^ — 0.65), the 
SHC keeps this minimum value unchanged. As shown 
in the following, this quantized SHC can be understood 
by the Z2-valued topological invariant associated with 
this quantum spin Hall phasc^^. When the Fermi en- 
ergy increases to touch the bottom of the middle band at 
e^=-0.65(=-l + A/2), then the SHC suddenly switches 
up and rapidly increases when the Fermi energy goes 
through the middle two bands. When ep increases at 
the top of the middle two bands (ei?=1.65), a^y arrives at 
the maximum value e/27r. Then as the Fermi energy ep 
continues to vary in the second gap region (shaded area 
1.65^eF^2), the SHC a^y keeps this maximum value 
unchanged. When the Fermi energy increases to touch 
the bottom of the upper band at ep=2, the SHC then 
suddenly switches down and rapidly decreases when the 
Fermi energy goes through the upper band. Finally the 
SHC a^y decreases to disappear when the three spin- 
degenerate bulk bands are all fully occupied. 

Topologically, the quantum spin Hall insulating state 
at 1/3 or 2/3 filling can be seen by calculating a selective 
Zo-valued invariant which is related to the parity 
eigenvalues ^2m(ri) of the 2m-th occupied energy band 
at the four time-reversal-invariant momenta F^. Very 
recently, Guo and Frana^^ have numerically calculated 
the eigenstate of Tir; and they found that three ^'s are 
positive and one is negative. Although which of the four 
^'s is negative depends on the choice of the inversion 
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FIG. 3: (Color online), (a) Sketch of the 2D kagome lattice 
strip with two edges along the y direction. The red, blue, and 
green dots are used to distinguish the independent sites A, 
B, and C, respectively, (b) Energy spectrum of the kagome 
lattice strip (A'i,=40) with the Hamiltonian 7i=7io + Ti-so- 
The spin-orbit coupling strength is set as Aso^O.lt. The red 
and blue lines represent the edge states localized at the down 
and up edges of the system, respectively. And the circle and 
triangle label the up and down spins, respectively, (c) The 
Riemann surface of the Bloch function corresponding to 1/3 
filling. The purple and orange curves correspond to spin-up 
and spin-down channels, respectively. 



center, the product Ili£^{Ti)—{—lY is independent of this 
choice and determines the nontrivial Z2 invariant 1/=!. 
That confirms the kagome lattice system to be a quantum 
spin Hall TBI at 1/3 (or 2/3) filhng. 

On the other hand, the topological aspect of the quan- 
tized spin Hall phase can be distinguished by the differ- 
ence between the winding numbers of the spin-up and 
spin-down edge states across the holes of the complex- 
energy Riemann surface, Is=I^ — ^j.— . The SHC is then 
given by aly—Ig {e/iir). Using this topological index Is 
we have investigated the quantum spin Hall effect in the 
Kane-Mele graphene model^". For the present kagome 
lattice model we can also study the quantum spin Hall 
effect in terms of this topological winding index. For 
this purpose, let us first numerically diagonalize the to- 
tal Hamiltonian Ti. using the strip geometry. For con- 
venience and without loss of generality, we suppose the 
system has two edges in the y direction while keeping 
infinite in the x direction [see Fig. EJJa)]. The num- 
ber of sites A (or B,C) in the y direction is chosen to 
be Ny=AO. The calculated energy spectrum is drawn in 
Fig. [S^b). From this figure one can clearly see that there 
are spin edge states occurring in each energy gap. These 
gapless edge states in the truncated kagome lattice are 
topologically stable against random-potential perturba- 
tion, provided that the perturbation is small compared 
to the bulk gaps. We have numerically confirmed this 
fact. The Riemann surface of Bloch function is plotted 
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FIG. 4: (Color online). The SHC as a function of the Fermi 
energy tp with different values of the Rashba coefflcient. The 
black, red, blue, and pink lines correspond to Aij=0, 0.02t, 
O.lf, and 0.2t, respectively. The intrinsic SOI strength is set 
as Aso^O.lt. 

in Fig. Eljc) for 1/3 filling. According to Ref;^, the 
winding number of spin-up (spin-down) edge state in the 
lower gap is I^=—l (Jj^ = l), which gives /s=— 2 at 1/3 
filling. That means the SHC in this phase is quantized 
as a'^y=—{e/2-K). The Riemann surface of Bloch function 
at 2/3 filling are the same as that at 1/3 filling, except 
that the directions of the curves corresponding to differ- 
ent spin channels are inverse. So the winding number of 
spin-up (spin-down) edge state in the upper gap is /|=1 
(/j^=— 1), which gives /s=2. The corresponding SHC at 
2/3 filhng is then quantized as (J^y=e/2'K. These conclu- 
sions are consistent with those calculated by using the 
Kubo formula ([7]) [see Fig. (He)]. Note that although 
the non-trivial Z2 invariant v=l in the bulk analysis can 
confirm the quantum spin Hall TBI phase, it does not 
provide the information on the sign of the SHC. In con- 
trast, our winding- number analysis can resolve this sign 
at different fillings (i.e., (Tjj^=^(e/27r) at 1/3 and 2/3 
filling, respectively). 

The topological properties of the system in the insulat- 
ing state are stable even when the Rashba SOI is consid- 
ered. To clearly see this fact, we numerically calculated 
the SHC a^y with the Kubo formula ([7]) when the total 
Hamiltonian includes the Rashba SOI^^ 

^ cL^(s ^ d,j);^Cj<,2, (8) 
{y)cri(T2 

where Xr is the Rashba coefficient and is a vector 
along the bond the electron traverses going from site j 
to i. The calculated SHC is drawn in Fig. 4 for different 
values of A/?. Clearly, one can see that the amplitude 
of SHC in the insulating phase keeps e/27r unchanged 
when the Rashba SOI strength 0^Afl,/i<0.1. When the 
Rashba SOI is sufficiently large (for example \ii—Q.2t), 
the amplitude of SHC will depart little from the quan- 
tized value. However, the topology of this insulating state 
keeps unchanged, unless the bulk gaps disappear. 

In the following let us consider the spin-charge sep- 
aration in the 2D kagome lattice. If we reconnect the 
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FIG. 5: (Color online), (a) Two edges of the 2D kagome lat- 
tice are reconnected with weaker bonds (orange dotted lines) 
and its corresponding energy spectrum (b). (c) The restored 
bonds have a sign reversal along the right half bonds (violet 
dotted lines) and its corresponding energy spectrum (d). In 
these figures, the system size is set as Ny=AQ and the spin- 
orbit coupling Aso=0.1f. The bonds connected two edges are 
sot as 0.25 times of other ones. 



two edges of the kagome lattice shown in Fig. [3]^ a) with 
weaker bonds, two small gaps then reappear in the edge 
spectrum [see Figs. El^a) and (b)]. Topological excita- 
tions (edge solitons) are created by reversing the sign of 
the reconnected bonds along the right half row [see Fig. 
[5l[c)]. In this manner two defects with tt flux are intro- 
duced in the present kagome lattice. As a result, four 
degenerate in-gap spin states localized around these two 
defects are formed in each bulk gap [see Figs. [DJd)]. The 
corresponding energies of these in-gap states are e——t 
and 1.7t, respectively. Note that in the square lattice 
version of the Kane-Mele modeli^, the in-gap modes are 
precisely at zero energy, while in the present kagome lat- 
tice model the in-gap modes are no longer at zero energy. 

At 1/3 or 2/3 filling, occupation (unoccupation) of 
these in-gap modes leads to an excess (deficit) of 1/2 
fermion number per spin and per defectlS. Four differ- 
ent types of solitons with the following quantum num- 
bers are obtained when these in-gap modes are occupied 
by different ways: the chargeon /_|_(i/2)|,+(i/2)j, (charge 
-e, Sz=0); the holon /_(i/2)T,-(i/2)i (charge e, 5^=0); 
the two spinous /+(i/2)t ,-(1/2)1 (charge 0, Sz = ^) and 
/-(i/2)T,+(i/2)i (charge 0, 5*2=— |). Here the subscript -I- 
(— ) represents that the in-gap mode is filled (empty) and 
t (i) labels the spin mode. For example, /+(i/2)t,-(i/2)x 
means that the up-spin mode is filled and the down- 
spin mode is empty. For further illustration, the charge- 
density distribution for the chargeon /+(i/2)t,+(i/2)x is 
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FIG. 6: (Color online), (a) The charge density of a chargeon 
state and (b) the spin density of a spinon on a 24 x 24 lattice 
with periodic boundary condition. The chargon at coordinate 
(6, 6) is /+(i/2)T,+{i/2)i, it has charge — e and spin Sz=0. The 
spinon at coordinate (18, 18) is /+(i/2)t, -(1/2)1, it has charge 
and spin Sz~fi/2. 

calculated and shown in Fig. Efa), while the spin-density 
distribution for the spinon /_|_(i/2)t,-(i/2)x is plotted in 
[S]Jb). Here, we use a 24x24 lattice for calculation. 

The quantum statistics of these spin-charge sepa- 
rated solitons can be readily seen by using anyon fu- 
sion argumen t ^'^i^^ , which is based on the observation 
that the bound states of fractional excitation acquire 
non-trivial Berry phases on adiabatic exchange. Since 
the spin-up and spin-down bands in the present case 
decouple, so without loss of generality, let us consider 
a bound state of two identical spin-up solitons f (1/2)1 
(or /_(i/2)|), which carries charge — e (e) and flux 2tt^0 
thus is a fermion. According to the anyon fusion rule, 
which states that the exchange phase of a particle 
formed by combining n identical anyons with exchange 
phase is Q=n^9, one easily obtain that the exchange 
phase between two identical solitons /(r/2)T (/-(r/2)T) 
should be 1/4 that of fermions, i.e., ^(/(i/2)t, /(i/2)t) 
(or ^(/_(i/2)t, /-(i/2)t)=='=''''/4- Next let us consider a 
bound state of an /(i/2)t and /_(i/2)t soliton, which 
carries charge and should be a boson. Then the 



exchange phase between solitons /(i/2)t ^^id /_(i/2)t 
is given by 6'(/(i/2)t , /_(i/2)t)=T7i"/4. Since the spin- 
down band is the Hermitian conjugate of the spin-up 
band one obtain Oifati, fa2l)=-0{faih 1021)- Substi- 
tuting these results into the formula calculating the 
exchange phase between spin-charge separated solitons 

^(/aiTftii /a2T/32i)=^(/aiT' /a2T)+^(//3ii, /fti), o^c im- 
mediately concludes that the spinous, holon, and char- 
geon are all bosons. However, each spinon has nontrivial 
mutual exchange phase tt with the chargeon and holon. 

Before ending this paper, we would like to point out 
that the ways creating the spin-charge separated solitons 
in the kagome lattice are not the only one. Recently Guo 
and Franz^^ suggested one possible realization of frac- 
tional charges ±e/3 or ±2e/3 by trimerizing the kagome 
lattice, hke the ID way proposed by Su and Schrieffer^ 
Based on this possible realization, if we introduce the in- 
trinsic SOI in the considered Hamiltonian, then the spin- 
charge separated excitations are expected to appear with 
fractional charges ±e/3 and spins zLfi/2 and ±fi/6. We 
leave this and the other interesting related problems to 
future study. 

In summary, we have theoretically studied the quan- 
tum spin Hall effect and the spin-charge separation in 
the 2D kagome lattice. By using the topological winding 
numbers of the spin-edge states on the complex-energy 
Riemann surface, we have obtained that the SHC is quan- 
tized as ±(e/27r) when the system is in the insulating 
phase, which is consistent with the Z2 topological invari- 
ant analysis and the numerical linear-response calcula- 
tion. Furthermore, we have constructed the spin-charge 
separated solitons in kagome lattice by connecting the 
system's boundaries in twist. The quantum statistics of 
these solitons has also been discussed. 
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